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Abstract 

We calculate nuclear effect in Drell-Yan massive lepton-pair production in 
terms of parton multiple scattering in Quantum Chromodynamics (QCD). We 
present the nuclear modification to inclusive Drell-Yan cross section da/dQ"^ 
in terms of multiparton correlation functions. By extracting the size of the 
correlation functions from measured Drell-Yan transverse momentum broad- 
ening in nuclear medium, we determine the nuclear modification at 0{as/Q'^). 
We find that the nuclear modification strongly enhances the inclusive Drell- 
Yan cross section in the intermediate mass region (IMR): 1.5 < Q < 2.5 GeV. 
We argue that the Drell-Yan process is responsible for the majority of the 
observed dimuon enhancement in the IMR in the Pb-Pb collisions at CERN 
SPS. 

PACS Numbers: 24.85.+p, 12.38.-t, 13.85.Qk, 25.40.Ve 
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The inclusive Drell-Yan massive lepton-pair production in hadron-nucleus and nucleus- 
nucleus collisions is an excellent laboratory for theoretical and experimental investigations of 
the role of nuclear medium in strong interaction dynamics. The massive lepton-pair produc- 
tion is also a channel for studying nuclear medium effect in production of quarkonium states 
and intermediate vector bosons. Recent measurements of massive lepton-pair production in 
both hadron-nucleus and nucleus-nucleus collisions show a very little nuclear dependence in 
Drell-Yan continuum, da/dQ, for lepton pair mass Q larger than 4 GeV Because of 

this lack of nuclear dependence in Drell-Yan continuum, the reported strong ^/ip suppression 
by the NA50 Collaboration was normalized to the Drell-Yan continuum in the lepton pair 
mass range 2.9 < Q < 4.5 GeV 0. 

Recently, it was shown by the NA38 and NA50 Collaborations that muon pair production 
for dimuon invariant mass between the and the ^/ip m heavy nucleus-nucleus collisions 
exceeds the expectation based on a linear extrapolation of the p — A sources with the product 
of the mass numbers of the projectile and target nuclei |Q. The excess increases with the 
number of participant nucleons, and the ratio between the observed dimuon yield and the 
expected sources reaches a factor of 2 for central Pb-Pb collisions. There have been a lot of 
effort to attribute such an excess to the sources other than the Drell-Yan continuum, such 
as the enhancement of open charm production [^], thermal dimuons production |^, and 
secondary meson- meson scattering in nuclear medium 0. As shown in Ref. [Q, the Drell- 
Yan continuum is the dominant source of the dilepton production in this relevant mass 
range in Pb-Pb collisions with 158 GeV/nucleon beam. An enhancement in the Drell-Yan 
continuum is much more effective than all other sources for interpreting the observed excess. 
Since the ^/ip cross section in Pb-Pb collisions was normalized to the measured Drell-Yan 
continuum in the same experiment, it is important to understand the nuclear dependence 
of the Drell-Yan continuum in the same mass range [Q. It is the purpose of this Letter to 
calculate the nuclear modification to the Drell-Yan continuum in the mass range between 
the and the J/V', which is reorganized as the intermediate mass region (IMR). 

Since we are interested in the Drell-Yan continuum, or in general, Drell-Yan cross section 
at low invariant mass, Q, we ignore contributions from the intermediate vector boson Z. 
Therefore, the inclusive Drell-Yan cross section in a collision between hadrons, h{p)+h'{p') — »• 
7*(— > ll{Q)) + X, can be expressed in terms of the cross section for production of an 
unpolarized virtual photon of the same invariant mass 

where S = {p + p'Y is the cm. energy of the collision. The Wfj,i^{Q) in Eq. (|lD is a hadronic 
tensor that satisfies electromagnetic (EM) current conservation, Q^W'^'^^Q) = W^" {Q)Qi, = 
0. Because of the EM current conservation of VF^iy(Q); the polarization tensor for the 
virtual photon, P'^^(Q) = Ea=t,l ^T{Q)(^\{Q) in Eq. (|1|), is often replaced by -g''\ The 
T and L here represent the transverse and longitudinal polarization of the virtual photon, 
respectively. 

When energy exchange of the collision, Q, is much larger than the typical momentum 
scales of initial-state hadron wave functions, the Drell-Yan cross section, da/d'^Q, is dom- 
inated by a single hard collision between a quark and an antiquark from two incoming 
hadrons, as sketched in Fig. |l](a). Corrections due to interactions involving extra partons 
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from incoming hadrons are power suppressed by the ratio of momentum scale of the hadron 
wave functions and the energy exchange of the hard coUision. Therefore, leading power 
contribution to the Drell-Yan cross section can be factorized as [ITOll 



where superscript (s) indicates a single hard scattering and J2a,b run over all parton flavors. 
The 0's are parton distributions and the factorization scale dependence are suppressed. The 
daab/d'^Q represents the perturbatively calculable rate of partonic collision. At the lowest 
order of quark- ant iquark annihilation, the inclusive Drell-Yan cross section is given by [jll| 



^ j dy (t)q/h{x) (t)qih'{x') (3) 



dQ2 y<^sQ 
where the parton momentum fractions are given by 

x = ;^e+^ and ^' = ^e-^ (4) 

with rapidity y. The high order corrections in powers of to the inclusive cross section 
in Eq. are well-defined in QCD perturbation theory |12]. For hadron-hadron collisions. 



Drell-Yan cross sections calculated in perturbative QCD are consistent with all existing data 



for Q > 4 GeV [13 



In hadron-nucleus and nucleus-nucleus collisions, more partons are available at a given 
impact parameter. Consider the Drell-Yan process in hadron-nucleus collisions. Before 
the hard collision of producing the lepton-pair, partons from different nucleons can either 
interact between themselves, as sketched in Fig. |l](b), or interact with the incoming parton, 
as sketched in Fig. |l|(c). After taking a square of the scattering amplitude in Fig. |l|(b), the 
leading contribution to the Drell-Yan cross section comes from the phase space where the 
parton momentum k is pinched, which corresponds to the situation in which the interactions 
between partons from different nucleons took place "long" before the short-distance hard 
collision. That is, only a single parton from the nucleus participates the short- distance 
hard collision. The parton-parton interactions between nucleons are internal to the nucleus, 
and should not affect the dynamics of short- distance hard collision. Therefore, leading 
contributions from the multiple scattering process in Fig. [I|(b) to the Drell-Yan cross section 
should have the same factorized form as that in Eq. (Q), except the parton distributions of a 
free hadron are replaced by the parton distributions of a nucleus. Those interactions internal 
to the nucleus in Fig. |l](b) are responsible for the phenomena of nuclear shadowing and 
EMC effect, which represent the differences between parton distributions of a free nucleon 
and those inside a nucleus. The nuclear shadowing and EMC effect provide small nuclear 
size dependence to the inclusive Drell-Yan cross section, or Drell-Yan continuum, da/dQ. 

In this Letter, we concentrate on the contributions due to the type of parton multiple 
scattering shown in Fig. |l|(c). Unlike the parton momentum k in Fig. |l|(b), the parton mo- 
mentum k' in Fig. |I](c) is not pinched after we take a square of the scattering amplitude 



Therefore, the type of double scattering in Fig. |l](c) is a part of the hard collision for pro- 
ducing the lepton-pair. Because of the extra initial-state interaction, physical contributions 
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from the double scattering are suppressed by a power of Q^, known as the high twist con- 
tributions [0. Although they are normally small in hadron-hadron collisions, such power 
suppressed contributions can be important in hadron-nucleus and nucleus-nucleus collisions, 
because of enhancement of nuclear size [|l^ . 

Before presenting our calculation of the nuclear modification to the Drell-Yan continuum 
caused by the type of multiple scattering in Fig. ^c), we estimate the significance of such 
power suppressed modification. Let the inclusive Drell-Yan cross section in the collisions 
between nucleus A and nucleus B be approximated as 



dQ^ dQ^ dQ^ 

^AB^[1 + RMQ)], (5) 

where superscript (d) represents the double scattering - a source of leading power corrections. 
The da]fj^/dQ^ is the single scattering cross section given in Eq. (|^) with parton distributions 
of free hadron replaced by the effective parton distributions of a nucleon inside a nucleus. 
The Rab{Q) represents a ratio of the double scattering over the single scattering, and has 
the following dependence 

where (Lab) oc A^^^ + B^/^ is an effective medium length. Since da/dQ = 2Qdcr / dQ"^ , the 
Rab{Q) in Eq- (H) is also equal to the ratio of the Drell-Yan continuum. If we assume 
that the Rab{Q) is positive and is about 5% at Q = 4 GeV for the Drell-Yan continuum 
in hadron-nucleus collision, which is within the experimental uncertainty of dilepton yield 
at that energy 0], the ratio Rab{Q) will yield almost 50% enhancement to the Drell-Yan 
continuum in the central region of Pb-Pb collisions aX Q = 2 GeV. Such an enhancement 
could explain the majority of the dimuon excess observed in Pb-Pb collisions. 

In order to quantitatively understand the strong nuclear enhancement to the Drell-Yan 
continuum in the IMR, we have to address the following two questions: (1) can the leading 
power corrections to the Drell-Yan cross section be systematically calculated in QCD per- 
turbation theory, and (2) is the sign of the nuclear modification to the Drell-Yan continuum 
positive. As we will show below, answers to these two questions are both positive. 



As a direct consequence of the generalized factorization theorem , the leading power 



corrections to the inclusive Drell-Yan cross section or Drell-Yan continuum can be systemat- 



ically calculated in QCD perturbation theory |15]. For the following derivation, we choose a 



frame in which the hadron (or nucleus) A is moving along the +2;-axis, and the B is moving 
in —z direction. We introduce the average nucleon momenta, p = Pa/ A and p' = Pb/B. 
The ± components of any vector k are defined in terms of two lightlike vectors n and n are 

= k ■ n and k~ = k ■ n. These lightlike vectors are chosen such that n ■ n = 1. In our 
choice of frame, the momentum p is up to a mass term proportional to n. 

At the leading order in Og, the leading power corrections to the Drell-Yan cross section 
are given by the Feynman diagrams in Fig. |^, where the quark lines with a short bar are 
special quark propagators, which are equal to the normal quark propagators with the leading 
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twist contributions removed |T^. According to the generahzed factorization theorem, the 
leading power corrections to the Drell-Yan cross section in hadron-nucleus collisions can be 



factorized as ITS 



w!,t\Q) = E / dx'(f)g/H'{x') f dxdxidx2 

q '' 

X TgA{x, xi, X2) H^^^J{x, xi, X2; x'; Q) , (7) 

where runs over all quark and antiquark flavors. In the light-cone gauge of strong 
interaction, gluon field strength, F^°'{y~) = n'^dpA°'{y~), and the hadronic matrix element 
TqAix,Xi,X2) is given by ^ 

dy- dy^ rfy^ 



TqAix,Xi,X2) 



1 



1 



,x — Xi/ \x — X2/ J 2tt 271 27r 



(8) 



The partonic part H^^J in Eq. (|^) is given by the Feynman diagrams in Fig. ^ in which 
quark (antiquark) lines are contracted by ^7 ■ p (|7 ■ p') and gluon indices are contracted 
with \dap. The transverse polarization tensor dap = —ga/3 + na^p + Uafip. 

Because of EM current conservation, we can calculate the partonic hard part by con- 
tracting the H^^J with either —g'^^ or the photon polarization tensor P^"{Q) in any gauge 
choice of EM interaction. Although the final answer is the same, contributions from different 
diagrams to the final answer are different. In addition, with an explicit polarization tensor, 
we can keep track the contributions to the different polarization states of the virtual photon. 
In the following, we present a result calculated in the light-cone gauge of EM interaction. 
In this gauge, the virtual photon's polarization tensor is given by 



With 



= 2PfiQ)- 

Q+fifj. _^ Q-j^fj. in Q^j- frame, we have 
1 



n 



Q 

Pl\Q) ■ 



(Q) 



d^" , and 



{Q) 



2Q- 



(9) 



(10) 



2"" ' g^ 

Contracting the partonic part Hj^j) with the polarization tensors in Eq. (|10|), we find that 
only one diagram in Fig. |^(a) contributes to the longitudinally polarized virtual photon 
state, and two diagrams in Figs, ^(e) and (f) contribute to the transversely polarized virtual 
photon state. We obtain 



HS{Q)Pr{Q) = el 



9' 



XX 



7 



X S{Q+ - xp+) 6{Q- - x'p'-) 6\Qt) 



X X 



1 

5{Q 



XX 



7 (2^)'^ 



X2P^ 



XiP'^ 



X2 - Xi 



Xi - X2 



6{Q- -x'p'-)6\Q 



T) 



(11) 



(12) 
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where g is the couphng constant of strong interaction and (1/18) is the overall color factor. 
Total partonic part to the inclusive Drell-Yan cross section is equal to the longitudinal 
contribution in Eq. (|lT|) plus twice of the transverse contribution in Eq. (0). If we contract 
the HjfJ with 
different diagrams. 

As pointed out in Ref. 



gfj-i^^ we get the exact same total contributions, except they come from 







the perturbatively calculated hard part in the generalized 
factorization formula is independent of the structure - in particular the size - of the nuclear 
target. We need to find the medium size enhancement due to multiple scattering from the 
matrix element in Eq. (^. The y~ integrals in Eq. (|^) parameterize the distance between 
the gluon and quark. Generally, the integrals over the distance cannot grow with the 
size of the target because of oscillations of the exponential factors, unless the integrations 
of parton momentum fractions lead to vanishing exponential factors. This can be achieved 



we 



if the integrations of parton momentum fractions are dominated by unpinched poles ||T6 
Combine the partonic parts in Eqs. (|TT]) and (|T^ and the matrix element in Eq. (|^) 
find that the dxi and dx2 integrations are dominated by the region where Xi ~ x and X2 ~ x 
from the poles. That the poles are in the matrix element TgA and do not have explicit ie 
is due to our choice of light-cone gauge for the strong interaction. Within the light-cone 
gauge, we can either assign the ie according the prescription introduced in the Appendix 



of Ref. [|T9| or carefully keep track of the ie in the partonic part. For example, we left the 
following factor 

X — Xi \ / X — X2 



Xi + le/ \x — X2 — le 

from the right-hand-side of Eq. ([TTl). The numerator cancels the artificial poles in the matrix 
element and leaves the poles with correct ie to the partonic hard part. We can also work 
in a covariant gauge of strong interaction without worrying about the artificial poles caused 

and the same answer was derived. 
■ Xi + ie) and 1/ (x — X2 — ie) to carry out the dxi and dx2 



by the choice of light-cone gauge [0, 

Using the unpinched poles 1/ (x - 
integrations in Eq. (|^), and substituting the hadronic tensor W'^ into Eq. (0), we obtain 
the leading order medium size enhanced power corrections to the inclusive Drell-Yan cross 
section in hadron-nucleus collisions. 



da 



id) 
Ah' 



dQ^ 




2 T 



qF 



X 



(t^g/h'ix) 



(13) 



with X and x' given in Eq. (^). The term proportional to Tgp (its derivative) contributes 
to the longitudinally (transversely) polarized virtual photon state. The twist-4 quark-gluon 
correlation function Tgp^x) is defined as [pO| , pl[1 



Tgpix) 



dy 

27r 



dyi dy2 
2-K 



e{-y2)d{y- 



X (P^|V^(0)^F+"(y2-)F„+(i/r)V^(y-)|PA). 



(14) 



In deriving Eq. (|T^), we used the following approximation in the limit Xi — > x and X2 —>■ x, 

1 1 



S{X2 - Q+/P+) ^ 5(xi - Q+/P+) 



X2 - Xi 



Xi - X2 



X — Xi + ie X — X2 — ie 
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^ 6'{x - ( ^— ^ -) + - (15) 

where represents the terms that have less than two poles. 

At the leading order in a^/Q^, the medium enhanced power corrections to the inclusive 
Drell-Yan cross section in nucleus-nucleus collisions are given by 

^^(rf) A„^'^^ A^i^) 

^^^^ 

with the nucleon- nucleus contributions given in Eq. (p!3|). 

In order to know the sign and magnitude of the medium enhanced power corrections, 
we need to know both the correlation function, Tqp{x), and its derivative. Although the 
correlation function is nonperturbative, the generalized factorization theorem requires it 
to be a universal function. Therefore, we can find the information of Tqp{x) from other 
observables. 

In Ref. , Drell-Yan lepton-pair's transverse momentum broadening was calculated in 
terms of parton multiple scattering in perturbative QCD. It was found that the broadening is 
proportional to the same quark-gluon correlation Tqp{x). Using data on Drell-Yan transverse 



momentum broadening in hadron-nucleus collisions, and a model for Tgp{x) from Ref. 

Tqp{x) ^ A^/' (f)q/r^{x) , (17) 



it was found |^ that ~ 0.01 GeV^. The same model was also used to evaluate 2-jet 



momentum imbalance in photon-nucleus collision [2C], and a larger value of (~ 0.05 GeV^) 
was extracted from Fermilab data [^]. Differences between these two observables were 



discussed in Ref. iT6 



Taking the same model for the quark-gluon correlation function in Eq. ([T7|), and substi- 
tuting Eq. (|T6|) into Eq. (|]), we obtain 



E.elJdy {A^^ [2 0,/;v(x) -a:^0,/jv(x)] (l)gfr^ix') + B'/^[x ^ x']} 



Eq / dy (j)q/N{x) (j)q/N{x') 



(18) 



where runs over all quark and antiquark flavors, and x and x' are given in Eq. (^). 

To be consistent with our leading order calculation in a^, we use leading order parton 
distributions CTEQ5L for the numerical evaluation of the ratio Rab{Q) HHIl- Although the 
value of the A^ is not well fixed, we will use A^ = 0.01 GeV^ to calculate the Rab{Q), because 
it was extracted from a similar process - Drell-Yan transverse momentum broadening. The 
uncertainty in the value of A^ corresponds to a shift in overall normalization of the power 
corrections, we present the ratio Rab{Q) as a function of dilepton's invariant mass Q in 
Fig. 0. The solid line is for Pb-Pb collision at \fS = 17.8 GeV/nucleon with cm. rapidity 
coverage < ycm < 1, while the dashed line is for p-W collision at y/S = 30 GeV/nucleon 
with cm. rapidity coverage —0.52 < ycm < 0.48. 

From the dashed line in Fig. ^, the nuclear enhancement to the Drell-Yan continuum in 
the p — W collision can be as large as 20% at Q ~ 2 GeV. This is not in any contradiction 
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with the data on the opposite-sign dilepton production from NA38 and NA50 collaborations, 
which show a good agreement with the expectation from all sources . As shown in Ref . , 
a majority of lepton pairs in this IMR in hadron-nucleus collisions at = 30 GeV come 
from open charm DD decay. Actually, a simultaneous fit of the 4 set p — A data shows that 
the ratio of opposite-sign dileptons from the open charm decay and those directly from the 
Drell-Yan process is as large as 4.2 ± 0.9 at = 30 GeV If. Therefore, a 20% increase of 
lepton pairs from the nuclear enhancement of the Drell-Yan continuum corresponds to less 
than 5% uncertainty in lepton pairs from DD decay, which is within the error bar of the 
fitting. 

On the other hand, as shown in Fig. 9 of Ref. , the expected opposite-sign lepton pairs 
from the DD decay are less than the pairs from the Drell-Yan continuum in the Pb-Pb 
collision. This is because total production rate of open charm DD pairs are much more 
suppressed than that of the Drell-Yan continuum when the collision energy decreases from 
= 30 GeV for p - A to = 17.8 GeV for the Pb-Pb collision. With the relative 
dominance of the Drell-Yan continuum in this IMR, a 50% enhancement of the Drell-Yan 
continuum at Q ~ 2 GeV from the solid line in Fig. ^ is a very significant change for the 
expected dilepton sources. Because of the mass threshold for producing open charm DD 
pairs, we expect the power corrections to the dilepton source from open charm decay to be 
much smaller. 

We also emphasize that when the dilepton's invariant mass Q is much less than 2 GeV, 
the leading order power corrections, as shown in Fig. |^, become so large that higher order 
power corrections have to be considered for a reliable QCD prediction. In addition, there 
should be high order corrections in as to the leading power contributions as well as the 
power corrections. Since we are interested in the relative size of these two contributions, we 
expect the ratio Rab{Q) to be relatively stable when we add high order corrections in a^- 
Finally, we point out that sign of the power corrections is not necessary to be positive as 
what we find here. Actually, as shown in Ref. the power corrections to the longitudinal 



and transverse structure functions in deeply inelastic scattering carry an opposite sign. 

In conclusion, we calculated the nuclear size enhanced power corrections to the inclusive 
Drell-Yan cross section in the IMR in hadron-nucleus and nucleus-nucleus collisions. We 
find that the medium size enhanced power corrections significantly enhance the dilepton 
production in the IMR, in particular, in the Pb-Pb collision. Therefore, in order to under- 
stand the measured nuclear enhancement of lepton pairs in the IMR at CERN SPS, a new 
fit for the expected sources including the power corrections to the Drell-Yan continuum will 
be very important. 

This work was supported in part by the U.S. Department of Energy under Grant Nos. 
DE-FG02-86ER40251 and DE-FG02-87ER40731. 
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FIGURES 




FIG. 1. Sketch for Drell-Yan process in hadron-hadron and hadron-nucleus collisions: (a) single 
hard scattering, (b) multiple interaction internal to the nucleus, (c) multiple interaction affects the 
hard scattering. 




(c) (d) (f) 




(g) (h) 



FIG. 2. Lowest order Feynman diagrams contribute to the leading power corrections to 
Drell-Yan cross section. 
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FIG. 3. The ratio Rab{Q) in Eq. (^) as a function of dilepton mass Q in the Pb-Pb (solid) and 
p — W (dashed) cohisions. 
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